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Abstract 

Equations of motion for a real self-gravitating scalar field in the background of a black hole with 
negative cosmological constant were solved numerically. We obtain a sequence of static axisym- 
metric solutions representing thick domain wall cosmological black hole systems, depending on the 
mass of black hole, cosmological parameter and the parameter binding black hole mass with the 
width of the domain wall. For the case of extremal cosmological black hole the expulsion of scalar 
field from the black hole strongly depends on it. 
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I. INTRODUCTION 



The early Universe and the phase transitions there provide signs of the high-energy phe- 

n 

nomena which are beyond the range of contemporary accelerators [l]. Such topological 
defects as cosmic strings attracted a great interests and were widely studied in literature. 
Black hole cosmic string configurations studies revealed the evidences that cosmic string 
could thread it or could be expelled from black hole. In the case of extreme Reissner- 
Nordstrom (RN) black holes it was found that there was a range of the black hole param- 
ete. fo.. „h,cK *e expuls.o. of the vo.ex .ooU place fl B flfl B Q 
theory being the low-energy limit of the string theory it was also justified 
a vortex could be treated by the remote observer as a hair on the black hole. Moreover, 
extreme dilaton black hole always expel the Higgs field, one has to do with the so-called 



In dilaton gravity 
^ ^ 'lO, 11| that 



Meissner effect. The problem of vortices in de Sitter background was analyzed in Ref. |12j |. 
while the behaviour of Abelian Higgs vortex solutions in Schwarzschild anti de Sitter (AdS), 
Kerr, Kerr-AdS and Reissner- Nordstrom AdS (RN-AdS) was studied in Refs. jl^l ll^. The 
vortex solution for Abelian Hi ggs field equations in the background of four-dimensional black 
string was considered in Ref. [15 1. 

The idea that the Universe is embedded in higher-dimensional spacetime acquires much 
attention. The resurgence is motivated by the possibility of resolving the hierarchy prob- 



lem 



17[, i.e., the difference in magnitudes of the Planck scale and the electroweak scales. 



Also the most promising candidate for a unified theory of Nature, superstring theory pre- 
dicts the existence of the so-called D-branes which renew in turn the idea of hrane worlds. 
From this theory point of view the model of brane world stems from the idea presented 
by Horava and Witten 



m 



[19^ . Namely, the strong coupling limit of x heterotic 
string theory at low energy is described by eleven-dimensional supergravity with eleventh 
dimension compactified on an orbifold with Z2 symmetry. The two boundaries of spacetime 
are ten-dimensional planes, to which gauge theories are confined. Next 0, 21, 22 1 it was 
argued that six of the eleven dimensions could be consistently compactified and in the limit 
spacetime looked five-dimensional with four-dimensional boundary brane. 

The studies of the interplay between black holes and domain walls (branes) acquire 
more attention. The problem of stability of a Nambu-Goto membrane in RN-dS space- 
time was studied in j22|. The gravitationally interacting system of a thick domain wall 
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and Schwarschild black hole was considered in j2J, |25[ . Emparan et al. |2a| elaborated the 
problem of a black hole on a topological domain wall. 

In j^^] the dilaton black hole-domain wall system was studied analytically and it was 
revealed that for the extremal dilaton black hole one had to do with the expulsion of the 
scalar field from the black hole (the so-called Meissner effect). The numerical studies of 



domain wall in the spacetime of dilaton black hole was considered in Ref. |28[, where the 
thickness of the domain wall and a potential of the scalar field ip^ and sine-Gordon were 
taken into account. In the case of a real self- interacting scalar field in the background of 
RN solution it was found that in the extremal case there was a parameter depending on 
black hole mass and the width of the domain wall which constituted the upper limit for 
the expulsion to occur [29]. Dynamics of domain walls intersecting black holes, taking into 
account the evolution of the system during the separation was elaborated in 30|. 

The importance of introducing the cosmological into considerations comes not only from 
the theoretical point of view but also for the observational results of our Universe. Small 
cosmological constant may be an alternative to the dark energy in explanation the observed a 
type la suprenova acceleration. On the other hand, as far as a negative cosmological constant 
is concerned a tremendous interest has focused on issues related to the AdS spacetime. One 
of them is the AdS / CFT (conformal field theory) correspondence which states that conformal 
field theories in d-dimension are described in terms of supergravity or string theory on the 
product spacetime consisting with asjTiiptotically AdSd+i and compact manifold, providing 
that there are relations between data on the boundary Rd of the AdSd+i and the data in 
the bulk AdSd+i 

Our paper will be devoted to the behaviour of the AdS black hole thick domain wall sys- 
tem. A domain wall will be simulated by a self-interacting scalar field. We shall numerically 
analyze the behaviour of scalar fields for various value of cosmological parameter. The brief 
outline of our paper is the following. The next section is devoted to the basic equations of 
the considered problem. In Sec. IllllllVl we presented the boundary conditions of the problem 
as well as the numerical analysis of the equations of motion for the two cases of potentials 
with discrete sets of minima. Namely, we shall analyze the ip^ and sine-Gordon potentials. 
The conclusions and discussions will be presented in Sec. El 
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II. THE BASIC EQUATIONS OF THE PROBLEM 



In our paper we shall consider a spherically symmetric static black hole with negative 
cosmological constant. The metric of which is written as follows: 

= - (i - ^ + ^ - r^) A= + /. _ a % _ A,.i + ^^('*'' + '"^^-^^ 

where we define A = ^ and Q is the charge of the considered black hole. In what follows we 
shall denote by V{r) = 1 — ^ + ^ — 4r^. The condition V{r) = is a quadratic algebraic 
equation for r. As was mentioned in the standard closed-form for the above mentioned 
four roots is rather lengthy and not especially illuminating. In what follows the AdS black 
hole spacetime will be the background metric on which one solves the domain wall Eqs. of 
motion. The domain wall will be simulated by a self-interacting scalar field. 

A general matter Lagrangian with real Higgs field and the symmetry breaking potential 
of the form will be taken into account 

/:dn. = -^v^^v^^y^-uiy^), (2) 

where by (p we have denoted the real Higgs field, while by U{(f) the symmetry breaking 
potential having a discrete set of degenerate minima. The energy-momentum tensor for the 
domain wall yields 

Tijiv) = -)j^gij^m^^"^^ - U {ip)gij + ViipVjip. (3) 

For the convenience we scale out parameters via transformation X = ip/rj and e = SttGi]'^. 
The parameter e represents the gravitational strength and is connected with the gravitational 
interaction of the Higgs field. Defining V{X) = where Vp = Xrj'^ we arrive at the 

following expression: 



SnGCdw = ^ 



(4) 



where w = y s-kGVf '^^P^^s^^ts the inverse mass of the scalar after symmetry breaking, which 
also characterize the width of the wall defect within the theory under consideration. Having 
in mind (jlj) the equations for X field may be written as follows: 
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In our considerations we take into account for two cases of potentials with a discrete set of 
degenerate minima, namely 

the ip^ potential described by the following equation: 

A 



(6) 



and the sine-Gordon potential of the form 



U2{^) = \v'[l + cosiifi/r])]. 



(7) 



III. THE BOUNDARY CONDITIONS 



In the background of the black hole spacetime the equation of motion for the scalar field 
X implies 



y - 2Mr + Q2 - ^r^)drX 



+ 



1 



sin 9 



do 



sin 9deX 



1 dV 
w^dX' 



(8) 



Consequently, having in mind relation Q one can define the following quantity for the scalar 
field (p: 



E 



V r 1 



,2/^_ 2M^ g _ A 

r 3 



- V{X). 



(9) 



As we have to do with non-asymptotically flat spacetime it cannot be interpreted as an 
energy density of scalar field. On the horizon of the black hole one has the following boundary 
conditions: 



V\rBH)drX 



1 



r=rBH~ 



r\jj sin 9 



de 



sin 9dQX 



+ 



1 dV 
'^dX' 



(10) 



As in Refs. |22l] we shall restrict our investigations to the case when the core of the 
domain wall is located in the equatorial plane, i.e., 9 = 7t/2. One ought also to impose the 
Dirichlet boundary conditions at the equatorial plane as follows: 



X 



0, 



(11) 



and the regularity condition for scalar field on the symmetry axis requiring the Neumann 
boundary condition on z-axis in the form as 

dX , 



09 



0. 



(12) 
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In order to find the solution of tlie problem our first task will be to establish the asymptoti- 
cally behaviour of the scalar field at the large distances. Further it will be necessary to find 
the numerical solution of its equation of motion. 

In this case at large distances from the horizon we demand that the black hole domain 
wall solution ought to be solution of the domain wall in AdS or dS spacetime. Defining the 
metric as follows: 

ds^ = -f 1 + ^^dt" + ^^^^ + r\de^ + sin^ ^V), (13) 

{} + P ) 

where = — 1 we have for AdS spacetime, on the other hand A; = 1 is responsible for dS 
one, the equation of motion for the scalar field X yields 

(l + '^)d^;X + H(l + ?|:!)S,.A- + is^Y + loots c)„A- - = 0. (14) 



As in Ref. [Ij] we set p = r cos 9 and the scalar field X is a function of p. By virtue of this 
equations of motion may be written in the form as 

Let us solve the equation for the potential V{X) = \{X^ — 1)^. The analitycal solution of 
Eq. is not known, but we shall search for the solution for magnitude of scalar field X{p) 
at constant large distances, i.e., for p —>■ oo. It follows in particular that X{p oo) = 1 = 
Xq. Equation (fT3j) is approximately satisfied by Xq ~ 1, where Xq is the minimum of the 
potential. Denoting fluctuations about this minimum by ip{p) we have the following: 

X(p)=Xo + ^(p). (16) 

It can be verified by expanding the scalar field X in the equation of motion that it reduces 
to the form 

kp^^^pp + AkptP^p + — ^ = 0. (17) 

In deriving Eq. (fTTj) one has neglected terms of order of unity in coefficients in first and 
second terms with respect to the terms involving derivatives of scalar field connected with 



£1 
P ■ 



-3+V9-4/3 
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Thus, Eq. fll7|) for large p —>■ oo has the approximate solution in the form of 

X(p)~Xo[l- 
6 



(18) 



where (3 = I'^/w^k, in our case P < 0. The same procedure apphed to the sine-Gordon 
potential V{X) = 1 + cosX case reveals the following relation: 

-2+^/4^-, 



1 - 



(19) 



IV. NUMERICAL CALCULATIONS 



A. Cosmological AdS black hole — (j)'^ potential 



Having laid out for our considerations, we proceed to describe the numerical analysis of 
the problem. First, we convert to equation of motion (jHJ to the form more suitable for finite 
differencing, namely one has the following: 

dr [fr' - 2Mr + + ^rA drX] -^de {smOdeX) - (x^ - l) = 0, (20) 

or making the differentiation with respect to r coordinate lead us to 

(r'' - 2Mr + + ^r^) d^X + 2 (^r - M + ^r^^ drX (21) 

+deeX + -^deX - ^X (x^ - l) = 0. 
tant^ ^ ' 

Next, we define the quantities 

1 . 

h = r 



a = r''- 2Mr + + -r^ 



M+-r\ 



(22) 



0. 



(23) 



It enables us to rewrite Eq. (j22j) in the form as follows: 

adrrX + 2bdrX + dggX + -dgX -X (X^ - 1 

tan 6* ^ ' 

Using the notation from our previous work |28,J Eq. (|^^ can be rewritten in the finite 
differences form. It implies 

Y 9ir„„-u]r„ Y ir„ 

(24) 



a — — 77! V 2b- 



(Ar) 



2Xoo + Xo 



1 Xo+ — Xo 
tan^ 2Ae 



w 



2Ar 

;Xoo (Xqo — = 0, 



or it can be rewritten in a more suitable form for our considerations. Consequently, Eq. (|25|1 
yields the result that: 



a b 
(Ar)2 Ar 



X+0 + 



+ 



2A^tan^ 



Xo- + 



(Ar)2 Ar 
2a 



X-0 + 



+ 



1 



2A9tan9 



Xo^ 



(25) 



(Ar)2 (A^)^ 



'{Xqq - I] 



X 



00 



0. 



We solve the equation ()26p using simultaneously over-relaxation (SOR) method 0. A 
general form of the second-order elliptic equation, finite differenced, yields 



y4x+o + Bx^o + Cxo+ + Dxq^ + Ex, 



00 



(26) 



Equation pHjl can be brought to the form of relation ()26|) with the following coefficients for 
the SOR method for our problem: 



E 



C 
D 



A 



B 



a b 
a b 



(Ar)2 Ar' 
1 1 



(Ae)2 2Aetan^fc^ 
1 1 



(A^)2 2A0tan^fc 



(Ar)2 (A^)2 ti;2 



F = 0. 



(27) 
(28) 
(29) 
(30) 

(31) 
(32) 



1. Boundary condition on the event horizon of black hole 



On the event horizon of the black hole under consideration a = and it can be verified 
that equation of motion (jHJ takes the form 



2 ( TEH- M+-r%H ]drX,r=r,, 



1 

tan 6 



0, 



(33) 



which yields the finite difference equation for our problem 

Oh ^+o--^oo , Xo+-2Xoo + Xo- 1 Xo+ - Xo- ^"eh ^ ( v'i A n /'Q/^^ 
2^^^ Ar + {Kef + tan^ 2A^^ " " = ^^^^ 

It can seen that by virtue of relation ()34j) one obtains the following: 



2b 



EH 



Ar 



X+0 + 



2A^tan^ 
2a 



^0+ + 

2 



1 



1 



(Ae)2 2A^tan^ 



EH I v2 



(Ar)2 (A^)^ 



{X. 



00 



Xo- + 
Xqo = 0, 



(35) 



where we have denoted by bEH the following expression: 

2 o 

bEH = Teh - M + -^r^jj. 



(36) 
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Now, the coefficients for the SOR method follows: 



^0 - ^, (37) 

Bo = 0, (38) 

^° = JKey + 2M^e ^^^^ 
_ 1 1 

i^o = 0. (42) 



2. Extremal black hole 

For the extremal black hole both a = and 6 = on the event horizon of the black hole. 
The field equation on the horizon decouples from the rest of the space and takes the form 

deeX + -^deX - aX (x^ - l) = 0, (43) 
tant^ ^ ^ 

where a = r^^/w"^. 



As was shown in Ref. j29| for a < 2 we have a field expulsion from the extremal RN 
black hole and one can set X = as a boundary condition on the event horizon. For a > 2 
equation ()43|1 must be solved prior to the calculation for the rest of the grid. For the case 
of the extremal black hole we shall adopt two point boundary relaxation method 



B. Cosmological AdS black hole — Sine-Gordon potential 

For the sine-Gordon potential the only difference in the SOR method are the E and F 
coefficients, which in this case take the form 

2a 2 

2 

F = -^sinXoo. (45) 

On the other hand, on the event horizon we have the following coefficients: 

2hEH 2 

2 

Fo = -^sinXoo. (47) 
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C. Anti-Nariai black holes 



Now we pay attention to the anti-Nariai black holes. The metric can be obtained by 
Ginsparg-Perry procedure 3^ from near-extreme AdS black hole. The AdS black hole that 
generates anti-Nariai one is of hyperbolic type. The anti-Nariai solution can be get when 
the the black hole horizon approaches the cosmological horizon. Its line element may be 
written in the form as js^] 

ds^ = -N{R)dT^ + + nl{d9^ + sinh ^'#'), (48) 



where N{R) = -1 + and 

1 



l + i^o 2 ^0-1^2 



A = = -^K- (49) 



For anti-Nariai case one has to replace the definitions (j22|l and (j36|) with 

a = RlN{R), b = RlN'{R), (50) 

and 

bEH = RIN'{Reh). (51) 
Then the SOR coefficient are the same as in the previous case with only one exception of 

^ = " " ~ 5^^°° ~ ^° = ~ " 5^^°° ~ ^^^^ 

for (p^ potential, and 

F = Fo = -^sinXoo, (53) 

for sine-Gordon potential. To conclude this section, we comment briefly on results of our 
numerical investigations. First we have studied the behaviour of the black hole domain wall 
system due to changes of cosmological parameter. Next we pay special attention to extreme 
black holes and domain walls. These systems are of special interests because of the possibility 
of expulsion of the scalar field from the black hole (the so-called Meissner effetct). Namely, 
Figure n represents the results of numerical integrations of equations of motion for the AdS 
black hole with 0^-potential and the mass M = 1.0, whereas the cosmological parameter 
/ = 1. The width of the domain wall was taken to be = 1. We depicted the values of the 
scalar field X and the quantity E = In Figs. llHUwe take into account the growth of 
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FIG. 1: The field X (left panels) and the quantity E (right panels) for the 4>'^ potential and the 
cosmological AdS black hole. Isolines on bottom panels are drawn for 0.2, 0.4, 0.6 and 0.8 for the 
field X and for —0.1, —0.2, —0.3 and —0.4 for the energy. Inlets in bottom plots show the value 
of the fields on the black hole horizon. Black hole has M = 1.0, Q = 0.1, / = 1.0 and the domain 
width is w = 1.0. 




FIG. 2: Same as Fig. d but for I = 2.0 

the cosmological parameter. Namely, in Fig. El we have / = 2, while in Fig. 01 one has / = 5 
and in Fig. jHZ = 10. The remaining parameters are the same as in Fig. ^ 

Figures OHl depict the solutions of equations of motion for scalar field X and show changes 
of E, for AdS black hole with sine-Gordon potential. The mass of the black hole under con- 
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FIG. 3: Same as Fig. HJ but for / = 5.0 




FIG. 4: Same as Fig. d but for I = 10.0 

sideration one has equal to M = 1.0 and the cosmological parameters are chosen respectively 
to be 1, 2, 5 and 10. The width of the domain wall was taken to be 1. On the other hand, 
Figures lUHTUl show the behaviour of the scalar field X on the black hole axis both for 0"^ and 
sine-Gordon potentials. 

Fi gure ITTl present s values of the ratio Q^/M^ for which the black hole under considerations 
becomes extreme. We consider it as a function of the cosmological parameter / and plots 
are made for different values of the black hole masses, i.e., M = 1, 10, 100, 1000. 

In Figs. IT^llTKl we plotted the values of X and the parameter E for extremal black holes. 
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FIG. 5: The field X (left panels) and the quantity E (right panels) for the sine-Gordon potential 
and the cosmological AdS black hole. Isolines on bottom panels are drawn for 0.27r, OAtt, O.Gvr and 
O.Svr for the field X and for —0.5, —1.5, —2.5 and —3.5 for the energy. Inlets in bottom plots show 
the value of the fields on the black hole horizon. Black hole has M = 1.0, Q = 0.0, / = 1.0 and the 
domain width is w = 1.0. 




0,40.81, 21, ( 




FIG. 6: Same as Fig. El but for I = 2.0 



It revealed that the field expulsion takes place only when the value of a parameter is smaller 
than 2.0. The same situation we ancountered in the case of RN black holes 2^. On the 
contrary, in the case of extremal dilaton black hole domain wall system one gets always 



expulsion of the scalar field 



271 |28|. It was revealed that the simplest generalization of 
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FIG. 7: Same as Fig. El but for / = 5.0 




FIG. 8: Same as Fig. but for / = 10.0 

Einstein-Maxwell theory by adding a massless dilaton field dramatically changes the struc- 
ture and properties of extremal black holes. Figure IT^ depicts the value of the scalar field X 
and E ioT a ^ 0.6 and for 0'^-potential. In Fig. we change the value of a. Now it is equal 
to 6.6. The results depicted in Figs. ITHITHl are valid for sine-Gordon potential and a = 0.6 
and 6.6, respectively. 

In Figs. [TMT7I we show the X and E values for anti-Nariai black hole with Rq = 1 and 
Kq = 1/2 both for 0^ and sine-Gordon potentials. In Ref. js^ the charged Nariai and 
anti-Nariai domain wall systems were studied. Due to the complication of equations of 
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0.98 - 



0.96 - 



0.94 - 



0.92 - 




2 4 6 8 10 

y 

FIG. 9: Values of the field X on the black hole axis for I = 1,2, 3, 5, 10 - solid curves from left to 
right, respectively. Dashed line represents the solution for I = inf. Rest of the parameters are the 
same as in Figure ^ 




2 4 6 8 10 

y 

FIG. 10: Same as Fig. IHl but for sine-Gordon potential. 

motion the simple arguments concerning the problem of expulsion of scalar field were given. 
However, the present analysis reveals the fact that there is no expulsion for the anti-Nariai 
case. Nevertheless, in the vicinity of the event horizon one obtains very small value of the 
field X. 
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0.5 1 1.5 2 2.5 3 3.5 4 
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FIG. 11: Q"^ /M"^ as a function of the cosmological parameter / for extreme black hole. The cm'ves 
are for M = 1, 10, 100, 1000, from left to right. 





r 



FIG. 12: The field X (left panels) and E (right panels) for the (f)^ potential and the extremal 
cosmological AdS black hole. Isolines on bottom panels are drawn for 0.2, 0.4, 0.6 and 0.8 for 
the field X and for —0.5, —1.5, —2.5 and —3.5 for E. Inlets in bottom plots show the value of the 
fields on the black hole horizon. Black hole has M = 1.0, the cosmological parameter / = 1.0 and 
a « 0.6. 

V. CONCLUSIONS 

Our considerations were devoted to studies of AdS black hole domain wall system. As 



in our previous works [2; 



29| the domain wall's equations of motion were built from a 
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FIG. 13: Same as Fig. [H but for a ^ 6.6. 




FIG. 14: Same as Fig. 1121 but for sine-Gordon potential. 

self-interacting real scalar field with a symmetry breaking potential having a discrete set 
of degenerate minima. In our considerations we took into account two kinds of potentials, 
i.e., 0^ and sine-Gordon potential. We used in numerical analysis simultaneously over- 
relaxation method 35| to solve equations of motion for scalar field X fulfilling the certain 
conditions on the black hole event horizon. The solutions of equations of motion depended 
on the parameter w = 1/ VXr] which was responsible for the domain wall thickness and the 
cosmological parameter /. We paid a special attention to extreme black hole solutions. As 
in the case of RN domain wall systems studied in Ref. j3| the expulsion of scalar fields 
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FIG. 15: Same as Fig. d but for a ^ 6.6. 




FIG. 16: X field configuration for anti-Nariai black hole Rq = 1.0 and Kq = 1/2 for (p^ potential. 
Left panel shows the value of the field, and right panel presents 0.2,0.4,0.6,0.8 contours. 




FIG. 17: Same as Fig. 1161 but for sine-Gordon potential. Contours on the right panel are for 
0.27r, O.lvr, 0.67r, O.Svr. 



18 



from AdS black hole is dependent on the parameter a = r^jj/w"^. The expulsion takes 
place when a < 2, while for a > 2 the domain wall can penetrate the extremal black hole. 
We confront the numerical studies of anti-Nariai black hole domain wall system to our very 
simple analytical analysis of the problem. Contrary to simple analytical arguments given 



in Ref. 



we revealed that numerical studies contradicts this suggestion. There is no 



expulsion of scalar field in this case. 
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